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Laterally coupled quantum dot molecules are studied using exact diagonalization techniques. We 
examine the two-electron singlet-triplet energy difference as a function of magnetic field strength 
and investigate the magnetization and vortex formation of two- and four-minima lateral quantum 
dot molecules. Special attention is paid to the analysis of how the distorted symmetry affects the 
properties of quantum-dot molecules. 
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I. INTRODUCTION 



The crossover from two-dimensional electron systems 
(2DES) to meso- and nanoscale quantum dots (QDs) is 
an interesting subject. In the infinite quantum Hall sys- 
tems the actual arrangement of impurities or disorder 
does not play a role, even if the presence of disorder- 
induced localized states is vital for the Hall plateaus to 
occuri In the few-electron QDs, the type of disorder is 
certainly an important issue. In the past, the majority of 
studies have concentrated on highly symmetric parabolic 
QDs without disorder. The rich spectrum of crossing 
energy levels as a function of magnetic field and strong 
interaction effects are nowadays rather well character- 
ized for the cases with a symmetric confinement poten- 
tial^ Recently, the focus has turned to understanding 
properties of QDs in less symmetric confinement poten- 
tials. For example, the rather simple far-infrared excita- 
tion spectrum of a purely parabolic QD is nowadays well 
understood^ whereas the lowered symmetry introduces 
new features in the spectrum whose interpretation is not 
straightforward at all.i&iiiLSiiLiSiil Moreover, the low- 
ered symmetry also gives rise to modified ground state 
properties such as level anticrossings and altered spin- 
phase diagram as a function of magnetic fieldi 1 ^ 3 " 

After Loss and DiVinccnzo proposition, 14 coupled 
quantum dots have gained interest due to possible re- 
alization as spin-qubit based quantum gates in quantum 
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In addition to coherent single-spin op- 



erations, the two-spin operations are sufficient for assem- 
bling any quantum computation. Recent experiments 
^ , have shown a remarkable success in characterizing the 
few-electron eigenlevelsfi 8 ^ 1 ^* 2 ^ approximating relaxation 
and time-averaged coherence times and mechanisms £L22k 
and reading single-spin or two-spin states2^24 of the QDs 
whereas the coherent manipulation of spin systems re- 
mained out of reach until very recent measurements on 
two-spin rotations^ 

In this paper we concentrate on two-electron quan- 
tum dot molecules. These molecules consist of laterally, 
closely coupled quantum dots. We treat correlation ef- 
fects between the electrons properly by directly diago- 
nalizing the Hamiltonian matrix in the many-body ba- 
sis (exact diagonalization technique). This allows direct 



access to the ground state energy levels and all excited 
states for both spin-singlet and spin-triplet states. We 
study these levels as well as singlet-triplet splitting and 
magnetizations as a function of magnetic field and dot- 
dot separation. We also analyze the properties of many- 
body wave functions in detail. 

The magnetic field dependence of the ground state en- 
ergy and singlet-triplet splitting in non-parabolic QDs 
have attracted recent interesti 12 i 13 i 18 i 26 i 27 i 28 i 29 i 30 i 31 i 32 i 33 
Magnetizations in QDs have been measured indirectly 
with transport measurements^ and recently with a di- 
rect technique with improved sensitivity^ 5 . For both mea- 
surements, semi-classical approaches cannot explain the 
results. Also the magnetizations of nanoscale QDs do not 
show non-equilibrium currents and de Haas-van Alphcn 
oscillations which are observed in 2DES and mesoscopic 
QDsi2& In the nanoscale QDs the quantum confine- 
ment and Coulomb interactions modify the system com- 
pared to the 2DES>S Theoretically, magnetization (at 
zero temperature) is straightforward to calculate as the 
derivative of the total energy with respect to magnetic 
field. Magnetizations have been calculated for a small 
number of electrons in a parabolic QD, in a square 
dot with a repulsive impurity^ as well as for anisotropic 
QDs. 2 ' and for self-assembled QDs and quantum rings. 39 
The magnetization curves have been calculated using 
density-functional theory for rectangular QDsSi and us- 
ing the Hartree approximation for other types of non- 
circular QDs. 40 A tight-binding model for 10-100 elec- 
trons in a single or two coupled QDs has been used to 
calculate magnetization curves. 41 

Calculations of vortices in QDs have also attracted 
much interest lately^Si 4 ^ 4 * 4 ^^ Even if the vortices are 
not directly experimentally observable, they reveal in- 
teresting properties of electron-electron correlations and 
of the structure of the wave function. The nucleation 
of vortices in QD systems could perhaps be observed by 
measuring magnetizations, where each peak would corre- 
spond to one vortex added in the system. However, the 
magnetization is difficult to measure for a small number 
of electrons, especially with direct methods. Moreover, 
in a non-circular symmetry, as in quantum dot molecules, 
and at high magnetic field strengths the magnetization 
curves may become more complicated. 

In our previous studies we have examined the proper- 
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ties of two-electron, two-minima quantum-dot molecules 
(QDM) in a magnetic field. The ground state as a 
function of magnetic field was found to have a highly 
non-trivial spin-phase diagram and a composite-particle 
structure of the wave function^ Also the calculated 
far-infrared absorption spectra in two-minima QDMs 
revealed clear deviations from the Kohn modes of a 
parabolic QD. Surprisingly, the interactions of the elec- 
trons smoothened the deviations instead of enhancing 
themi£ In Ref. H3 we briefly discuss some of the results 
of square-symmetric four-minima QDM. 

In this paper, we study in detail the properties of dif- 
ferent QDMs. Three different QDM confinements are 
studied thoroughly and their properties are compared 
to parabolic-confinement single QDs. First, we calcu- 
late measurable quantities such as energy eigenstates, 
singlet-triplet splittings and magnetizations as a func- 
tion of magnetic field. Secondly, non-measurable quan- 
tities, such as conditional densities, vortices, total densi- 
ties, and the most probable positions are used to analyze 
the nature of the interacting electrons in quantum states 
and also to analyze and understand the properties of the 
measurable quantities. This paper is an extension to 
our previous calculations of QDMspi 2 ^! We study a two- 
minima QDM (double dot), a square-symmetric four- 
minima QDM and a rectangular-symmetric four-minima 
QDM. The aim of this paper is to study how the con- 
finement potential affects the properties of interacting 
electrons in a low-symmetry QD. 

This paper is organized as follows: In Section[n]we ex- 
plain how the quantum-dot molecules are modeled and 
what kind of basis we use in the exact diagonalization 
method. We also discuss calculation of magnetizations 
and the conditional single-particle wave function which 
we use to locate the vortices and study the conditional 
density. In the following four sections we present our 
results. In Section IIHI we discuss properties of a sin- 
gle parabolic quantum dot, and in Sec. IIVI we analyze 
the properties of a double dot. In Sec. [3 we present re- 
sults for the square-symmetric four-minima quantum-dot 
molecule, and finally in Sec. the results for rectangu- 
lar four-minima quantum-dot molecule. The analysis of 
the results is presented in Section I V 111 The summary is 
given in Section lYlIII 



II. MODEL AND METHOD 



We model the two-electron QDM with the two- 
dimensional Hamiltonian 



H = 
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FIG. 1: Confinement potential of square-symmetric (L x 
L y — 5 nm) four-minima quantum dot molecule. 



where V c is the external confinement potential taken to 
be 



V c (r) 



-m 0Jn mm 



(2) 



where the coordinates are in two dimensions r = (x, y) 
and the L/s (Lj = (±L X7 ±L y j) give the positions of the 
minima of the QDM potential, and M is the number of 
minima. When Li = (0, 0) (and M = I) we have a single 
parabolic QD. With M = 2 and Li,2 = (±L X , 0) we get a 
double-dot potential. We also study four-minima QDM 
{AI = 4) with minima at four possibilities of (±L X , ±L y ) 
(see Fig. We study both square-symmetric [L x = 
Ly) and rectangular-symmetric (L x ^ L y ) four- minima 
QDMs. The confinement potential can also be written 
using the absolute values of x and y coordinates as 



V c (x,y) 
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For non-zero L x and L y , the perturbation to the 
parabolic potential comes from the linear terms of L x or 
L y containing also the absolute value of the associated 
coordinate. 

We use the GaAs material parameters m* /m e = 0.067 
and e = 12 .4, and the confinement strength huiQ = 3.0 
meV. This confinement corresponds to the harmonic os- 
cillator length of Iq = y/h/ui^m* w 20 nm. We concen- 
trate on closely coupled QDMs where L x<u < Iq. The 
magnetic field (in z direction) is included in the symmet- 
ric gauge by the vector potential A. The Hamiltonian 
of Eq. is spin-free, and the Zeeman energy can be 
included in the total energy afterwards (Ez — 9* H-bBSz 
with g* = —0.44 for GaAs). We disregard the threefold 
splitting of each triplet state [Sz = 0, ±1) and consider 
only the lowest energy one {Sz — !)• 

We drop the explicit spin-part of the wave function 
and expand the many-body wave function in symmetric 
functions for the spin-singlet state [S — 0) and anti- 
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FIG. 2: Relative error in the energy of parabolic two-electron 
QD at B = as a function of the basis size n x — n y from 3 
to 8, which corresponds to around 40-2000 many-body basis 
functions in the expansion. The relative angular momentum 
state of electrons is denoted with m. 

symmetric functions for the spin-triplet state (S = 1). 
*s(ri,r 2 ) =^a i ,j{<?!'i(r 1 )^(r2) 

+ (-l) s <Mr 2 )^-(r 1 )}, (4) 

where ctjj's are complex coefficients. The one-body basis 
functions 4>i(r) are 2D Gaussians. 



( r ) = x n *y ny e~ 
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(5) 



where n x and n y are positive integers. The complex co- 
efficient vector ai and the corresponding energy Ei are 
found from the generalized eigenvalue problem where the 
overlap and Hamiltonian matrix elements are calculated 
analytically. The matrix is diagonalized numerically. 

The basis is suitable for closely coupled QDs. At large 
distances and at high magnetic field we expect less accu- 
rate results. The accuracy may also depend on the sym- 
metry of the state. At zero magnetic field the parabolic 
two-electron QD can be modeled with a very good pre- 
cision by expanding the basis (in a given symmetry) in 
relative coordinates. In Fig. [21 we compare the energy 
of the very accurate solution and the one using our basis 
(for a parabolic QD) as a function of the basis size, where 
the maximum n x = n y ranges from 3 to 8. These values 
correspond to around 40 — 2000 many-body configura- 
tions in the expansion. States with m = 0, 1, 2, 3 refer to 
different relative angular momentum states. The relative 
error, even with the smallest basis studied n x = n y = 3, 
is less than one percent and decreases rapidly with the 
increasing basis size. The greatest error is found for the 
rn = state. 

The magnetization can be calculated as the derivative 
of the total energy with respect to magnetic field. It can 
be divided into to parts, paramagnetic and diamagnetic, 
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FIG. 3: One-body density (dotted line), two-body spin- 
singlet state (dashed line) and two-body spin-triplet state 
(solid line) along x axis. One of the electrons is fixed at 
the most probable position in the x axis (x*) and the condi- 
tional density is plotted for the other electron (|?/) c (r)| 2 = 
\<fsl(x,y),(x*,0)}\ 2 /\^ s [(-x*,0),(x*,0)}\ 2 ). The peaks on 
the left-hand side also indicates the most probable position, 
therefore by reflecting the peak position to the right-hand side 
of the x axis one can perceive the position of the fixed elec- 
tron, (a) and (b) show the densities of parabolic QD at two 
different magnetic field values (B — 1 and 8 T), (c) and (d) 
represent two-minima QDM with L x = 20 nm. The confine- 
ment potential, V c , is plotted with gray color on each figure. 



where the former is constant as a function of magnetic 
field, for a given angular momentum and spin state, and 
the latter depends linearly on magnetic field. The dia- 
magnetic contribution to the magnetization is also a mea- 
sure of the spatial extension of the ground stated 

Total electron density can be obtained by integrating 
one variable out from the two-body wave function 



n(ri) = J dr 2 |^ s (ri,r 2 )| 2 . 



(7) 



In practice we do not perform numerical integration. The 
density is directly calculated in our diagonalization code 
where the required matrix elements are calculated ana- 
lytically. 

We analyze the two-body wave function by construct- 
ing a conditional single-particle wave function 

*s[(x,y),(xt,i£)] 



Mr) = \Mr)\e l 



®s[(x$,yl), (£2.2/2)] 



(8) 



where one electron is fixed at position (xjjj , y|) and the 
density (IV'c(r)l 2 ) and phase {0 c (r)) of the other electron 
can be studied. One of the electrons is usually fixed at 
the most probable position (x^J/f), but we also analyze 
V> c (r) when the other electron is fixed at some other po- 
sition. The most probable positions of electrons (rj,^) 
are found by maximizing the absolute value of the wave 
function with respect to coordinates ri and T2'. 



max|*s(ri,r 2 )|^ 

ri ,r 2 



(9) 
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One should note that | C r ) 1 2 ^ s n °t normalized to one 
when integrated over the two-dimensional space because 
it describes the electron at position (x, y) on the con- 
dition that the other electron is fixed at (x^,^)- I n ~ 
stead, IV^r)! 2 is normalized so that it equals to one 
when x = x*,y = yj. Using the conditional single- 
particle wave function we can study the conditional den- 
sity ^(r)! 2 and the phase 9 c (r). 

To illustrate how the properties of the many-body wave 
function can be examined with the conditional single- 
particle wave function, we compare interacting two-body 
conditional densities to non-interacting two-body densi- 
ties in Fig. 13 The non-interacting two-body density is 
the same as the single-particle density (up to a normal- 
ization). We call it the one-body density hereafter. We 
plot the one-body, two-electron singlet {S = 0) and two- 
electron triplet (5 = 1) conditional single-particle densi- 
ties along x-axis. The other electron, in the two-electron 
systems, is fixed at the most probable position (x*) on 
the right-hand side of the x-axis. 

Fig. |3 (a) and (b) show conditional densities of the 
single parabolic QD at B = 1 and 8 T magnetic fields. 
The one-body density is located at the center since no 
correlation effects push it towards the edges of the dot. 
The peak of the triplet state is found further at the edge 
of the dot than the singlet peak since Pauli exclusion 
principle ensures that the electrons of the same spin are 
pushed further apart than the electrons with the opposite 
spins. Notice that the conditional density of the triplet 
state goes to zero where the other electron is fixed, just 
before x = 20 nm, whereas in the singlet state there is a 
finite probability to find the electron around the point of 
the fixed electron. 

Fig. 13(c) and (d) show the same data for L x — 20 nm 
two-minima QDM. In high magnetic fields and at large 
dot-dot separations the difference between singlet and 
triplet densities reduces. In QDMs, with a sufficiently 
large distance between the dots and in high magnetic field 
also the one-body density localizes into the individual 
dots (see Fig. 131(d)). 



III. PARABOLIC TWO-ELECTRON QUANTUM 
DOT (L = 0) 

We start our analysis from the single parabolic quan- 
tum dot. The two-electron parabolic QD is studied ex- 
tensively in the literature but presenting results here 
serves as a good starting point for understanding prop- 
erties of quantum dot molecules. 



A. Energy levels 

Energy levels of the parabolic QD are plotted in Fig.|3] 
as a function of magnetic field. Fig. 0] (a) shows non- 
interacting two-body energy levels, (b) the ten lowest 
levels for two-body spin-singlet states (S = 0), and (c) 



for two-body spin-triplet (S = 1) states. In (d) three 
lowest singlet and triplet levels are shown in the same 
plot. The non-interacting spectrum is obtained by oc- 
cupying two electrons in the Fock-Darwin energy lev- 
els. The first eigenvalue at zero field equals two times 
(N e = 2) the confinement potential (huio = 3 meV, 
E\{B = 0) = 3 + 3 = 6 meV) and the second level repre- 
sents one electron in the lowest Fock-Darwin level and the 
other electron in the next one (E2(B = 0) = 3 + 6 = 9 
meV). Many non-interacting energy levels are degener- 
ate, also as a function of magnetic field. (In a less sym- 
metric confinement, the degeneracies are lifted). Due to 
degeneracies, only six levels are seen in Fig. SI (a). If the 
interactions are included, the spectra become much more 
complicated and many more level crossings are observed. 
One can also see how the energy scale is modified. In 
the spin-singlet spectra the ground state energy is almost 
doubled if the Coulomb interaction is included. 

To see the crossing singlet and triplet states more 
clearly, we plot the three lowest energy levels of spin sin- 
glet (dashed line) and spin triplet (solid line) in Fig. 01(d) 
up to B — 6 T. In a weak magnetic field the ground state 
of the two-electron QD is spin-singlet (S = 0), which 
changes to spin-triplet (S — 1) as the magnetic field in- 
creases and then again to singlet and finally to triplet (at 
B « 6.3 T, not visible in Fig.0|(d)). 



B. Singlet-triplet splitting and magnetization 

In Fig. \5\ (a) we plot the energy difference of the lowest 
triplet and singlet states up to B = 8 T. Altering singlet 
and triplet states are also seen in higher magnetic fields 
with a decreasing energy difference between the states. 
However, if one includes the Zeeman term, the triplet 
state is favored over the singlet state at high B. Therefore 
the system becomes spin polarized. 

The transitions between the states can also be exam- 
ined from the magnetization curves plotted in Fig. [51(b). 
The non-interacting magnetization is a smooth curve as 
no crossings are present in the lowest energy level. The 
orbital angular momentum in the non-interacting two- 
electron ground state does not change, and thus only 
the diamagnetic effects are seen in the magnetization. 
The non-interacting electrons have a smaller spatial ex- 
tent of the wave function compared to the interacting 
electrons. Therefore, at low fields, when the response is 
purely diamagnetic, the magnetization curve of interact- 
ing electrons has a lower absolute value in Fig. |S] (b). 
The magnetization curve of interacting electrons shows 
abrupt increase of the otherwise smooth curve whenever 
two levels cross (see also Fig. SI (d)). The peaks in mag- 
netization are solely due to interactions. 
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FIG. 4: Ten lowest energy levels of parabolic QD (L x = 0, L y = 0) as a function of magnetic field of (a) non-interacting 
two-body, (b) two-body singlet state (S = 0) and (c) two-body triplet state (S = 1). Some of the states in (a) are degenerate, 
(d) three lowest energy levels of singlet (dashed line) and triplet (solid line) as a function of magnetic field up to B = 6 T for 
parabolic QD (L x = L y = 0). In (b) the first singlet ground state corresponds to angular momentum m = which changes to 
m — 2 at B ~ 2.7 T. In (c) the first triplet ground state is m = 1 and it changes to m = 3 at B ~ 5.8 T. In (d) the first ground 
state equals to m — singlet, then ground state is m = 1 triplet followed by m = 2 singlet and m = 3 triplet, where the latter 
is not visible as a ground state in (d). Zeeman energy is included in the triplet energies (Ez = —2 x 12.7B[T] ^eV in GaAs). 
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FIG. 5: Triplet-singlet energy difference in (a) and magnetiza- 
tion in (b) for parabolic QD (L x — L y — 0). The lower curve 
in (a) shows the singlet-triplet splitting with Zeeman energy 
included. The smooth curve in (b) represents the magneti- 
zation of two non-interacting electrons and the other curve 
shows the magnetization for two interacting electrons. At 
low magnetic field values, the ground state is the angular mo- 
mentum m — singlet, which shows as positive values in the 
triplet-singlet energy difference of (a) and as a smooth curve 
in (b) . When the first transition from m = singlet to m = 1 
triplet occurs triplet-singlet energy difference changes its sign 
from positive to negative and there appears a peak in magne- 
tization. Change of m — 1 triplet to m — 2 singlet and from 
m = 2 singlet to m = 3 triplet appear in the same way as 
peaks in magnetization and changes of sign in triplet-singlet 
splitting. Magnetization is given in the units of effective Bohr 
magnetons fj* B = eh/2m* (n* B = 0.87 meV/T for GaAs). 



C. Wave function analysis & vortices 

We will now analyze the two-body wave functions and 
study in more detail singlet-triplet transitions in the sin- 
gle QD. The first singlet-triplet transition can be be 
understood with the simple occupation of the lowest 
single-particle states: In the singlet state the two elec- 
trons occupy the lowest energy eigenstate with opposite 
spins (S — 0). As the magnetic field increases, the en- 
ergy difference between the lowest and the second low- 



est single-particle levels decreases. (See non-interacting 
energy levels in Fig. 21 ( a ))- At some point the ex- 
change energy in the spin-triplet state becomes larger 
than the energy difference between the adjacent energy 
levels. Thus, the singlet-triplet transition occurs and the 
adjacent eigenlevels are occupied with electrons of paral- 
lel spins (S = 1). 

However, the true solution of the two-electron QDM 
is much more complicated than the occupation of single- 
particle levels and inclusion of exchange energies. Inter- 
action between the electrons changes the situation dras- 
tically. This can already be seen by comparing the single- 
particle energy levels of Fig. (a) to singlet and triplet 
energy levels in (b) and (c). As a signature of complex 
many-body features, many singlet-triplet transitions are 
seen as a function of B. There are two trends competing 
in the ground state of a quantum dot when the magnetic 
field increases. The magnetic field squeezes the electron 
density towards the center of the dot and the Coulomb 
repulsion of electrons increases at the same time as the 
electron density is forced to a smaller volume. At some 
point it is favorable to change to a higher angular momen- 
tum ground state, which pushes electron density further 
apart and reduces the Coulomb energy. Therefore as a 
function of the magnetic field a series of different angular 
momentum states are seen. 

The altering singlet and triplet states can also be un- 
derstood in terms of composite particles of electrons and 
attached flux quanta^ The starting point for under- 
standing the ground state changes and flux quanta is to 
consider the two-electron parabolic QD (as discussed in 
this Section), which has an exact solution for the wave 
function of the form 

* = (xia + iyi2) m f(ri2)e-^ +r i )/2 , (10) 

where x\ 2 = X\—x%, y\i = y±— 2/2 and r i2 = |rj -r 2 | are 
the relative coordinates of the two electrons, m is the rela- 
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FIG. 6: (a)-(e) contours of conditional density \ip c (x, y)\ 2 and 
phase of the conditional wavefunction 8 c (x,y) in gray-scale 
for parabolic QD (L x = L y — 0). (White equals 9 C = and 
darkest gray 8 C = 2n). Magnetic field value and the spin type 
are plotted on top of each figure. The plus sign (+) indicates 
the position of the fixed electron and small circles indicate 
the positions of the vortices. In (f) contours of total electron 
density of the three-vortex triplet state are plotted in the 
background and positions of the vortices are solved when the 
fixed electron is in three different positions. Fixed electron is 
marked with the plus sign and vortices with circles. The most 
probable position is marked with a star, on the left-hand side 
for clarity. 



tive angular momentum and / is a correlation factor 
The zeros of the wave function (vortices in relative coor- 
dinates) are placed on Z12 = £12+22/12 = with a winding 
number given by the relative angular momentum to. In 
the first S = state the relative angular momentum of 
electrons is zero (to = 0). When the magnetic field in- 
creases the ground state changes to spin triplet 5 = 1, 
where the relative angular momentum of electrons equals 
one (to = 1) and in the second singlet state m = 2, and 
so on. With increasing magnetic field the relative angular 
momentum of electrons increases and altering singlet and 
triplet states are seen (if the Zeeman term is excluded). 
The transitions occur because in the states with large 
to the Coulomb repulsion becomes smaller at the cost of 
higher confinement and kinetic energies. As the increas- 
ing magnetic field squeezes electrons to a smaller area, it 
is favorable to move to larger to to minimize the total en- 
ergy. One should note that with even to the spatial part 
of the total wave function is symmetric and therefore the 
spin part should be antisymmetric (S — 0). With odd to 
the spin part is symmetric (S = 1). 

When the angular momentum increases, the increased 
rotation induces vortices in the system. As we have a 
many-body system, the rotation is a correlated motion of 
electrons and can be studied in the relative coordinates 
of electrons. Vortices can be found by locating the zeros 
of the wave function and studying the phase of the wave 
function when going around each of the zeros. As the 
vortices are seen in the relative coordinates, and are not 



visible in the density, we examine the conditional single- 
particle wave function ^> c (r) (of Eq. JSJ), where one elec- 
tron is fixed in the most probable position (on the x axis, 
the system is circular symmetric). The vortices are seen 
in the zeros of ip c . When the phase part, 9 Cl is integrated 
around a closed path encircling the zero, we obtain the 
winding number of the vortex (<f> 9 c (r)dr — to27t). 

In a parabolic two-electron QD vortices are automati- 
cally attached on top of the electrons, where the relative 
angular momentum to equals the winding number of a 
vortex. In a less symmetric potential the center of mass 
and relative variables do not decouple and one would ex- 
pect more complicated structures as can be seen in later 
sections. The simple form of the wave function in Eq. 
HI U| l is due to separation of the center of mass and rela- 
tive coordinates in the parabolic confinement. 

Calculated vortices and conditional densities of the sin- 
gle QD are shown in Fig. HJ1 The contours show the con- 
ditional electron density, \^p c \ 2 , and the gray-scale back- 
ground marks the phase of the conditional wave func- 
tion, 9 C , where the white color equals 9 C = and the 
darkest gray 9 C = 2ir. The positions of the vortices are 
marked with circles (o), and the other electron is fixed at 
the most probable position (r|) shown with a plus sign 
(+)-The lines of dark gray and white borders correspond 
to a sudden phase change of 27r if the line is crossed. The 
number of flux quanta attached to the electron (or the 
winding number of a vortex) can be determined by go- 
ing around the fixed electron position and calculating the 
total phase change (or counting the lines crossed in the 
figure) . 

In Fig. (a) the phase is constant (no vortices and 
no relative angular momentum) and the the probability 
density of the other electron is located on the left side 
because of the Coulomb repulsion. In (b) we find one 
vortex as one border of white and gray is crossed when 
the fixed electron is encircled. In (c) we find two vortices, 
in (d) three vortices and in (e) four vortices. Fig. H3 (a) 
corresponds to the first singlet with relative angular mo- 
mentum to = 0, (b) shows the first triplet with to = 1, 
(c) the second singlet with to = 2 and (d) the second 
triplet with to = 3. FigEl ( c ) would be the third singlet 
state but this is not a ground state if the Zeeman term 
is included. The conditional density shows how the elec- 
tron localizes to a smaller area when the magnetic field 
increases. We also notice the enhancement of interac- 
tion at high B where the density contours are contracted 
compared to the low-field conditional densities. 

The conditional density and phase are much more sen- 
sitive to the basis size than e.g. energy eigenvalues. In 
a parabolic QD the vortices should appear exactly on 
top of the fixed electron. However, in the two-, three- 
and four-vortex plots, in Fig. El the vortices are slightly 
displaced from the fixed electron position. The finite ba- 
sis expansion does not result in exactly correct vortex 
picture. On the other hand, the problem is not very se- 
rious because the error in energy is not large and vortex 
positions are not experimentally observable and we also 
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FIG. 7: Total electron density of the ground state at different 
magnetic field values for parabolic (L x = L y = 0) QD. In the 
left column are densities in 1/nm 2 and on the right column 
are contours of the densities. Dark regions in the contours 
mark high density and brighter regions low density. Notice 
that in the densities the scale in the x and y axes is changing 
and the height of the peak is increasing with magnetic field. 
The range of axes in the contour plots is kept fixed. The 
magnetic field value and the spin type of the ground state is 
plotted above each sub-figure. 



get correct winding numbers (angular momentum) if the 
pinned electron is encircled with a large enough radius. 

Vortex dynamics can be studied if we change the po- 
sition of the fixed electron. In Fig. (f) we show the 
positions of the vortices in three different places of the 
fixed electron for the three-vortex triplet state. The to- 
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FIG. 8: Most probable position in nm of singlet (S = 0) and 
triplet (S = 1) states for parabolic (L x = L v — nm) QDM. 
The magnetic field region where singlet is a ground state is 
marked with gray background color. 

tal electron density is shown in the background (see also 
Fig. [3(g) and (h)). In a parabolic QD the vortices should 
be on top of the fixed electron no matter where it is 
pinned. Our calculations result wrongly a small offset 
of vortices from the fixed electron. Even if the vortices 
are not exactly on top of the fixed electron, and actually 
have a greater offset at greater distances from the origin, 
they are seen to follow the fixed electron if its position 
is changed. This is a signature of the composite particle 
nature of electrons and vortices. 



D. Total electron density and the most probable 
position 

In Fig.[7|we plot ground state total electron densities of 
the parabolic QD. The magnetic field values are the same 
as in Fig. except that the four-vortex solution is not 
plotted since it is not a ground state. The range of both 
x and y axes and also the density in the z axis changes 
from (a) to (g) whereas in the contours the range of x and 
y is kept fixed. Dark regions in the contours mark high 
density and bright regions low density. As a function of 
magnetic field there forms a minimum in the center of 
the dot as the Coulomb repulsion forces electrons further 
apart. Note that even if there is a minimum in the total 
electron density, this is not a vortex. Vortices, in the 
composite particle picture, follow moving electrons and 
are seen in the relative coordinates of electrons. They 
describe the correlated motion of electrons and are not 
visible in an averaged-out quantity such as density. 

Another way to study the nature of the changing 
ground states is to plot the most probable position (see 
Eq. 0|) as a function of magnetic field. Fig. [H] shows 
the most probable position (r* ) for both spin-singlet and 
spin-triplet states as a function of magnetic field. Gray 
background color indicates the region of magnetic field 
where the singlet is a ground state and white background 
color corresponds to magnetic field regions where the spin 
triplet is a ground state. Each jump in the curves corre- 
sponds to a change in the angular momentum. When the 
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in the single QD, but the transition points change and 
the energy differences are smaller at greater distances 
between the dots than in the single QD. The transition 
points and regions of singlet and triplet states are plotted 
in Fig [5] (b) . We can also note that all transition points 
are shifted to lower B at large distances between the 
dots. If the Zeeman energy, that lowers the triplet energy, 
is included in the total energy the second singlet state 
disappears at greater L as can be seen in Fig. (c) and 
(d). The second singlet is only seen in a small region 
with very closely coupled QDs (L < 2.5 nm). Therefore, 
subsequent singlet states after the first one are not seen 
in the double dot if L > 2.5 nm. Similar results are 
seen in anisotropic QDs where the parabolic confinement 
of a single QD is elongated continuously to a wire-like 
confinement^ 



FIG. 9: Triplet-singlet energy difference (AE = E n - E n ) 
as a function of magnetic field in two-minima quantum dot 
molecule. The energy difference is plotted as a function of 
dot-dot separation and magnetic field in (a) without Zeeman 
energy and in (c) with the Zeeman energy included (AE = 
E^ + Ez — E^). In (b) and (d) the ground state regions of 
the singlet and triplet states are plotted as function of B and 
L, without and with Zeeman energy, respectively. 

singlet changes from m — to m = 2 state at B ps 2.7 
T the most probable position jumps to a higher value 
as well. See also energy level crossings in Fig. 0] (b). 
The outward relaxation, due to the increase of angular 
momentum, can be also seen in the density. Similar de- 
pendence is seen for the triplet state. First the most 
probable position decreases due to contracting electron 
density and then, at some point, it is favorable to move 
to a higher angular momentum state which relaxes the 
electron density outwards. In higher magnetic field we 
would see a sequence of transitions between increasing 
angular momentum states. 



IV. TWO-MINIMA QUANTUM DOT 
MOLECULE (L x 0) 

A. Singlet-triplet splitting as a function of L 

In this section we study two laterally coupled quantum 
dots. In a two-minima QDM, or double dot, we study the 
changes in the ground state spectrum when two QDs, on 
top of each other, are pulled apart laterally. In Fig. |U| 
(a) the energy difference of the lowest triplet and singlet 
states is plotted as a function of the inter-dot spacing 
and magnetic field. At L = we have a single parabolic 
QD and the curve coincides with Fig 0(a). When L^O 
we have a double dot. Let us now examine some general 
trends of the triplet-singlet energy difference as a function 
of dot-dot separation. At small magnetic field the ground 
state is a spin singlet, then triplet, and again singlet as 



B. Energy levels of L x = 5 nm double dot 

We choose one distance between the dots, L x = 5 nm, 
to study the properties of the double dot in more detail. 
We plot energy levels, singlet-triplet splitting, magne- 
tization, vortices, the most probable position and total 
electron density of this double dot in Figs. El - El The 
energy levels in Fig. 1101 are now modified, compared to 
the parabolic QD, due to the lower symmetry of the con- 
finement potential. The circular symmetry is no longer 
present. The lower symmetry shifts and splits degenerate 
levels. The non-interacting levels, in Fig. El( a )j split a t 
zero magnetic field and there is also a small anticrossing 
of levels, just barely visible in the figure. Also degenerate 
levels of the single QD (see Fig. 0] (a)) are now slightly 
displaced in energy, which is mostly seen as thicker lines 
in Fig. El (a). In the interacting spectra, (b) and (c), 
we see many anticrossings. Also the nature of the lowest 
level does not change abruptly with crossing levels as in 
the single QD, but instead we see anticrossing levels. For 
example, there are clear anticrossings in the double dot 
singlet states of Fig. El (b) whereas states cross in the 
singlet state of parabolic QD of Fig. 0] (b) . 

We plot three lowest singlet and triplet energy levels 
in the same figure to see the transition points and energy 
differences between the states more clearly (Fig. El (d)). 
The ground state is a singlet at small B, also in the dou- 
ble dot, and later it changes to triplet. The Zeeman term 
lowers the triplet energy enough so that no second singlet 
(ground) state is observed at higher B, even though the 
singlet energy becomes very close to the triplet energy at 
B w 5 T, as can be seen in Fig. El (d). The interesting 
anticrossing of spin singlet between B = 2 and 3 T is 
now an excited state as the triplet is the ground state. 
We also find anticrossing ground state levels in the spin 
triplet around B w 5.5 T, but the repulsion of levels is 
not so clear at high B. 
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FIG. 10: Energy levels of L x = 5, L y — two-minima QDM. See Fig.0]for details. 





FIG. 11: Triplet-singlet energy difference (a) and magnetiza- 
tion (b) for L x — 5, L y = two- minima QDM. See Fig. 0for 
details. 



Singlet-triplet splitting and magnetization of 
L x — 5 nm double dot 



The energy difference between the lowest triplet and 
singlet states as a function magnetic field in the L x = 5 
nm double dot is plotted in Fig. 1111 fa), and the magne- 
tization in Fig. 1111 (h). The sharp jump in the magneti- 
zation corresponds to the singlet-triplet transition. Even 
if the system is not circular symmetric, and angular mo- 
mentum is not a good quantum number, there is an in- 
crease of the expectation value of angular momentum at 
the transition. We can clearly see that the magnetiza- 
tion increases suddenly at the transition point. Around 
B s» 5.5 T there is a bump in the magnetization. This 
is exactly at the anticrossing point of the triplet state. 
Therefore the symmetry of the triplet state changes or 
the magnetic moments of the electrons change. This time 
it is not seen as an abrupt change but as a continuous 
one. Similar magnetization curves are seen in asymmet- 
ric QDs with a correct deformation from the parabolic 
confinement to a more wire-like confinement^ 



D. Vortices of L x — 5 nm double dot 

In the case of QDMs the states cannot be identified 
with angular momentum since it is not a good quantum 
number. However, we can still study vortices and condi- 
tional density of the double dot. We fix one electron at 




FIG. 12: (a)-(e) contours of conditional densities \ip c (x,y)\ 2 
and phase of the conditional wavefunction 8 c (x,y) in gray- 
scale for L x — 5, L y = two-minima QDM. See Fig.HJfor de- 
tails, (f ) contours of total electron density of the three- vortex 
triplet state are plotted in the background and positions of 
the vortices with the fixed electron in three different positions. 



the most probable position at r* = (x*, 0) and study the 
conditional single-particle wave function 



4>c(x,y) 



V s [(x,y),{x*,0)] 

0),(x*,0)]" 



(11) 



The most probable position of the two-minima QDM lies 
always on the x axis. In Fig.^](a), at low B, no vortices 
are found and the phase is constant. Contours are again 
localized to the left of the fixed electron having the max- 
imum at (— x*,0). The conditional density in a double 
dot is more localized compared to the conditional density 
of a single QD in Fig. © (a). The next plot, Fig. d (b) 
shows data for the triplet at B = 3.0 T with one vortex 
and (c) shows the singlet state at B = 5.7 T with two 
vortices, (d) the triplet at B — 7.5 T with three vortices, 
and (e) the singlet at B = 8.2 T with four vortices. The 
singlet state is not a ground state after the first singlet- 
triplet transition which means that two- and four-vortex 
solutions in (c) and (e) are only found as excited states. 
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FIG. 13: Most probable position in nm of singlet (S = 0) and 
triplet (S = 1) states for two-minima (L x = 5, L y = nm) 
QDM. Singlet ground state magnetic field region is marked 
with gray background color. 
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FIG. 14: Density of the ground state at different magnetic 
field values for two-minima (L x = 5, L y = nm) QDM. See 
Fig- for details. 



the vortices are seen in the vicinity of the fixed electron 
even if the conditional density is localized closer to the 
other dot and contour lines are rather circular though 
flattened in the x direction due to Coulomb repulsion. 
Despite the fact that the density of an electron is rather 
localized to one of the dots the two electrons move in a 
strongly correlated way in the double dot. The dynam- 
ics of vortices can be studied by changing the position of 
the fixed electron. In Fig. El (f) the electron is fixed in 
three different positions for the B = 7.5 T triplet ground 
state. The total electron density is plotted in the back- 
ground. The vortices follow the fixed electron which leads 
us to conclude that also in a non-symmetric potential 
flux quanta and electrons form composite particles. It 
is surprising to find composite particle solutions of elec- 
trons and flux quanta in the non-parabolic symmetry as 
well* 

The vortices are not exactly on top of the fixed elec- 
tron. This may be in part due to the finite-size basis 
expansion, but on the other hand there is no reason why 
they should be exactly on top of the fixed electron. In- 
teresting vortex clusters of a six-electron parabolic QD 
are studied in a recent article (Ref. l42f) and also in el- 
liptical and rectangular QDs in Ref. 45. In the double 
dot, the vortex patterns change continuously, where al- 
ways an extra pair of vortices approach the fixed electron 
from minus and plus infinity from the y axis lying in the 
same vertical line with the fixed electron as in Figs. El 
(c)-(e). 



E. Most probable position and total electron 
density of L x = 5 nm double dot 

The most probable position (r*) for the lowest singlet 
and triplet states is plotted in Fig. El Due to the anti- 
crossing levels also the singlet and triplet most probable 
positions change continuously. A similar outward relax- 
ation of the conditional density in the double dot is as- 
sociated to the approaching vortices as was seen for the 
parabolic QD associated with the sudden change in the 
angular momentum state. 

Ground state densities at three different magnetic field 
values are shown in Fig. 1141 They show how the electron 
density localizes into the two minima as the magnetic 
field increases. The interacting density shows two peaks 
also in the low-field regime whereas the non-interacting 
density still has just one maximum (see the dotted line 
in Fig. 13(c)). 



In the double dot vortices appear in the dot similarly 
as in the single QD. The vortices, on the other hand, are 
not exactly on top of the fixed electron, and they can be 
found at considerable distances from the fixed electron, 
as in Fig. El ( e )- Looking at the ground state at B = 7.5 
T in Fig. ll2l fd) in more detail, it is interesting to see that 



V. SQUARE-SYMMETRIC FOUR-MINIMA 
QUANTUM-DOT MOLECULE (L x = L y / 0) 

In this Section we study two electrons in lateral four- 
minima quantum-dot molecules. The minima are ar- 
ranged in the way that they form a square in the lateral 
direction. 
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FIG. 15: Triplet-singlet energy difference (AE = E n - E n ) 
as a function of magnetic field in square-symmetric (L = L x = 
L y ) four-minima quantum dot molecule. See Fig.[5]for details. 



A. Singlet-triplet splitting as a function of L 

We study singlet and triplet states as a function of 
magnetic field and dot-dot separation. Fig. El ( a ) shows 
altering singlet and triplet states as a function of mag- 
netic field. More frequent singlet-triplet changes are seen 
at greater separations between the dots (large L). Also 
notable is the large energy difference of the second sin- 
glet state to the triplet state. The second singlet also 
persists as a ground state to the greatest studied sep- 
aration L, which is not true in the double dot, if the 
Zeeman energy is included (Fig. El ( c ))- ln au studied 
separations L the magnetic field evolution of the Zeeman 
coupled four-minima QDM (Figs, (c) and (d)) follows 
the same pattern. At small magnetic field values the 
ground state is singlet, then triplet and again singlet in 
a small magnetic field window before the ground state 
changes to triplet permanently. However, with large sep- 
arations L the system becomes spin-polarized at lower 
magnetic field values, i.e. the border line of the second 
singlet and second triplet curves towards low-field region 
with increasing L. 

We will now analyze the rapid changes of the singlet 
and triplet states of four-minima QDM. Singlet-triplet 
(and triplet-singlet) transition points shift to lower mag- 
netic field values at greater L, where the area of the QDM 
(A) is effectively larger. If the transitions occur at effec- 
tively the same values of the magnetic flux (<f> = BA), 
the transitions should be seen at lower B when the area 
is larger. This explains why the border lines between sin- 
glet and triplet states curve towards lower B at greater 
separations between the dots. 

The second singlet is seen as a ground state in the 
Zeeman coupled system even at the very large distance 
of L = 20 nm where the perturbation from a purely 



parabolic potential is clear. However, in this type of 
square- or ring-like potential spatially symmetric states 
(singlet) are energetically more favorable than in elon- 
gated potentials (double dot) where in general the spa- 
tially antisymmetric states (triplet) are favored without 
paying too high price in Coulomb energy. Of course, as 
all energy scales are quite equal, the ground state is a del- 
icate balance between kinetic, confinement and Coulomb 
energies as a function of magnetic field. 



B. Energy levels of L x = L y — 5 nm QDM 

We will now focus on the L x = L y = 5 nm four-minima 
QDM. The energy levels of non-interacting two-body, sin- 
glet and triplet states are plotted in Fig. El ( a )-( c )> re- 
spectively. In the square-symmetric four-minima QDM 
the ground state levels do not anticross as in the dou- 
ble dot, instead crossing ground states of the same spin 
state are seen, as in the single QD. However, the four- 
minima QDM is not circular symmetric and anticrossings 
are still seen in the higher energy levels, which is not the 
case with circular symmetric parabolic QD. Notice that 
in the non-interacting spectra in Fig. 1161 (a) many lev- 
els are degenerate at B = for the four-minima QDM, 
whereas in the double dot the zero-field degeneracies are 
mostly lifted, see FigllUI (a). Yet, many levels that are 
degenerate in the single QD, Fig0] (a), as function of 
magnetic field, are slightly split in the non-interacting 
energy levels of four-minima QDM. 

The singlet and triplet energies can be seen in the same 
plot in Fig. 1161 (d) . The energy levels of the four- minima 
QDM in Fig 1 161 (d) look very similar to single QD energy 
levels of Fig|U(d). Even if the absolute values of energies 
and transition points arc different, the only notable dif- 
ferences between parabolic QD and four-minima QDM 
are the small bending at zero field of the third triplet 
level, near E w 12.9 meV, and small anticrossing of the 
uppermost triplet level at B w 1.3 T. 



C. Singlet-triplet splitting and magnetization of 

L x = L y = 5 nm QDM 

The energy difference between the lowest triplet and 
singlet and the magnetization are plotted in Fig. El 
Now, as no anticrossings of ground states are present, 
the triplet-singlet energy difference shows peaks. In the 
magnetization we also observe sharp peaks whenever two 
ground states cross. Similar results are seen in a square 
quantum dot with a repulsive impurity of Ref. |3^ where 
the magnetization of two electrons in the dot shows sharp 
transitions since no anticrossings in the ground states 
are present. On the other hand, in a square QD with 
eight electrons slightly rounded magnetization curves are 
seen^& Therefore, if the circular symmetry of the con- 
finement is broken, the magnetization depends on the 
symmetry of the confinement but also on the number of 
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FIG. 16: Energy levels of L x — L y — 5 nm four-minima QDM. See Fig. [I]for details. 
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FIG. 17: Triplet- singlet energy difference in (a) and magne- 



tization in (b) for L x 
details. 



5 QDM. See Fig. ^ for 
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FIG. 18: (a)-(e) Contours of conditional densities \ip c (x,y)\ 2 
and phase of the conditional wavefunction 8 c (x,y) in gray- 
scale for L x = 5, L y — 5 four- minima QDM. See Fig.HJfor de- 
tails, (f) contours of total electron density of the three- vortex 
triplet state are plotted in the background and positions of 
the vortices with the fixed electron in two different positions. 



the electrons in the QD. It may not be straightforward to 
draw any conclusion about the underlying potential from 
the magnetization curves. Anticrossings, on the other 
hand, are clear signatures of a broken circular symmetry. 
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FIG. 19: Most probable position in nm of singlet (S = 0) 
and triplet (S = 1) states for four-minima L x = L y = 5 nm 
QDM. Singlet ground state magnetic field region is marked 
with gray background color. 



D. Vortices of L x 



5 nm QDM 



We can identify the changes in the magnetization to 
the increasing number of vortices in the two-electron 
QDM or to the increase of the expectation value of rela- 
tive angular momentum of the electrons. Vortex patterns 
and conditional densities are shown in Fig. The most 
probable position is now found from the line connecting 
two minima diagonally. 



*s[(x,v),(x*,V*)] 
*s[(-x',-y*),(x*,y*)y 



y 



(12) 



At B — 0.7 T in (a) the conditional density is spread to 
the area of three unoccupied dots with a peak in the fur- 
thermost dot on the diagonal from the fixed electron. At 
high magnetic field the density becomes more localized 
closer to the most distant minimum, in the diagonal from 
the fixed electron. However, the contours show that the 
conditional density is not as circularly symmetric as in a 
double dot, but actually resembles more the conditional 
density of the single parabolic QD. The peak in the con- 
finement potential at the origin seems not to affect the 
conditional density considerably when compared to the 
single QD. At high B the Coulomb repulsion forces elec- 
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FIG. 20: Density of the ground state at different magnetic 
field values for four-minima L x = L y = 5 nm QDM. See 
Fig. |7|for details. 

tron density to the outer edges of the confinement, which 
might not result in very different results when compared 
to the single dot. However, the distance L x = L y = 5 
nm in the confinement is not particularly large and the 
perturbation from the parabolic confinement is not very 
large. On the other hand, altering singlet and triplet 
states persist to the greatest studied distance (L = 20 
nm) between the dots of L = 20 nm where the perturba- 
tion from the parabolic confinement is clear. 

The vortices appear in the four-dot QDM sequentially 
as a function of magnetic field. At low B the ground 
state is a singlet with no vortices, then a triplet with one 
vortex, a singlet with two vortices and then a triplet with 



three vortices. The singlet with four vortices in Fig. ITS1 
(c) is an excited state as the system becomes spin polar- 
ized after the two-vortex singlet state. The vortices are 
located in the diagonal going through the fixed electron 
position, see Figs. |^| (c)-(e). The vortices seem to be 
further away from the fixed electron (in the case of more 
than one vortex) than in the single QD. This is also true 
in the double dot. This could be identified to repulsion 
between the vortices but it is difficult to assess since the 
length scales are different due to different confinement 
strength and also the basis causes some errors. However, 
with six electrons in a parabolic confinement one can see 
a clear repulsion between the vortices 4£ 

Vortex dynamics of the three- vortex solution is studied 
by changing the position of the fixed electron in Fig. ^| 
(f). Total electron density of the same state is plotted 
in the background of Fig. ^] (f). Vortices are seen to 
follow the fixed electron also in the four-minima QDM. 
However, in the four-minima QDM the vortices are fur- 
ther away from the fixed electron as the distance from 
the origin increases. 



E. The most probable position and density of 
L x = L y = 5 nm QDM 

The most probable positions of the lowest singlet and 
triplet states of the four-minima QDM (Fig. I19f) show 
very similar behavior as the single QD. Only the most 
probable positions are on average roughly 2 nm greater 
at all field strengths compared to the single QD (size 
of the QDM is larger compared to single QD). Otherwise 
continuously decreasing r* shows a jump when the lowest 
singlet (or triplet) state changes. 

Ground state densities and contours are plotted in 
Fig. |20| Starting from a rather flat density at low fields, 
a hole begins to form in the center as the magnetic field 
is increasing. The electron density localizes into a nar- 
rowing ring around the origin. However, compared to the 
parabolic QD, there are peaks forming in the four corners 
of the density, instead of a smooth density ring. Also the 
density looks more square-like in all of the contours. 



VI. RECTANGULAR FOUR-MINIMA 
QUANTUM DOT MOLECULE (L x =£L V =£Q) 

In this Section we examine the triplet-singlet energy 
difference, energy eigenlevels, magnetization, vortices, 
the most probable positions, and densities of four-minima 
QDM with rectangular positioning of the QD minima in 
the lateral plane. 
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FIG. 21: Triplet-singlet energy difference (AE = E v ' - E rl ) 
as a function of magnetic field in rectangular-symmetric four- 
minima quantum dot molecule. L x is fixed to 5 nm and L y 
is varied from to 20 nm. Therefore at L y = we have 
L x = 5 nm double dot and at L y — 5 nm it is the square- 
symmetric four-minima QDM. The energy difference is plot- 
ted as a function of L y and magnetic field in (a) without 
Zeeman energy and in (c) with the Zeeman energy included 
(AE = E n + E z - E n ). In (b) and (d) the ground state 
regions of the singlet and triplet states are plotted as function 
of B and L y , without and with Zeeman energy, respectively. 



A. Singlet triplet splitting as function of L y for 
fixed L x — 5 nm rectangular QDM 

The triplet-singlet energy difference of the rectangular 
QDM is plotted in Fig. The distance between the 
minima in the x direction is fixed while the distance in the 
y direction is varied. We set L x = 5 nm and vary L y from 
zero to 20 nm. Therefore, at L y — we have a L x = 5 nm 
double dot, and at L y — 5 nm we have the L x = L y = 5 
nm square-symmetric four-minima QDM studied in the 
preceding section. The smooth surface in Fig. is due 
to anticrossing ground states similarly as in the double 
dot. The anticrossings in the lowest levels of the singlet 
and triplet states are again present if the symmetry is 
distorted from a square to a rectangular symmetry, see 
the energy levels in Fig. 1221 The only sharp peaks in 
Fig. 1211 (a) and (c) correspond to rectangular symmetric 
QDM at L x =L y = 5 nm. 

One can also see that as a function of magnetic field 
singlet and triplet states do not change as rapidly as in 
the square-symmetric four-minima QDM. Interestingly, 
the third singlet region terminates around L y w 15 nm. 
So at sufficiently large distance between the two L x = 5 
nm double dots the singlet state is no longer favorable 
even if the Zeeman term is excluded. In the case of a 
double dot in Ref. ^3 it was n °t possible to say whether 
the second singlet state would terminate at greater dis- 
tances between the dots, but for two double dots the third 



singlet region clearly terminates. 

If the Zeeman term is included (Figs. |^ (c) and (d)) 
the second singlet state can only be observed in a small 
region where the rectangular four-minima QDM is close 
to square symmetry (near L y = 5 nm). Actually the 
energy difference has its maximum, as function of L yi 
at the square symmetry. If we follow the energy differ- 
ence at zero magnetic field as a function of L y , it first 
increases reaching the maximum at L y = 5 nm and then 
it decreases again when L y is increased. 

The stability of the singlet states (and also triplet 
states) in the sgware-symmetric QDM can be understood 
from the relatively high energy of the triplet state (or 
singlet) near the peak in AE. In the square symme- 
try the degenerate energy levels at the crossing point are 
energetically very unfavorable. In rectangular symme- 
try degeneracies are lifted (anticrossings), which lowers 
the energy of the other spin type and also reduces the 
energy difference, AE. Thus the energy differences are 
always smaller in the rectangular symmetry when anti- 
crossings are present. The Jahn- Teller theorem states 
that any non-linear molecular system in a degenerate 
electronic state will be unstable and will undergo a dis- 
tortion to form a system of lower symmetry and lower 
energy, thereby removing the degeneracy 48 In a QDM, 
the system cannot, of course, lower the symmetry of the 
external confinement spontaneously to lift the degenera- 
cies, but the large triplet-singlet energy differences in 
the square-symmetric QDM can be understood via Jahn- 
Teller effect: If the symmetry is lowered, degeneracies are 
lifted and smaller triplet-singlet energy differences are ob- 
served. 



B. Energy levels of L x = 5, L y — 10 nm QDM 

We will now study rectangular L x — 5, L y — 10 nm 
QDM in more detail. Fig. |221 (b) an d (c) reveal many 
anticrossings in the interacting two-body spectra of L x = 
5, L y = 10 nm QDM, both in the ground states and also 
in the excited states. Many features look similar as in the 
double dot spectra in Fig. ^| but the anticrossing gaps 
are bigger here due to the greater separations between 
the dots (greater deviation from the circular symmetry). 
Fig. |221 (d) shows singlet and triplet energy levels in the 
same plot. The second singlet becomes very close to the 
triplet near B = 5T, but the triplet remains the ground 
state. 



C. Singlet-triplet splitting and magnetization of 

L x = 5, L y = 10 nm QDM 

The energy difference of triplet and singlet states is 
plotted in Fig. |23 (a). The magnetization in Fig. [2*^1 
(b) show first a sharp peak which corresponds to singlet 
triplet transition. The next change is from the one- vortex 
triplet to the three- vortex triplet and as these two states 
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FIG. 22: Energy levels of L x = 5, L y = 10 nm rectangular four-minima QDM. See Fig. ||| for details. 




FIG. 23: Triplet-singlet energy difference (a) and magneti- 
zation (b) in L x = 5, L y = 10 nm rectangular four-minima 
QDM. Magnetization is given in the units of effective Bohr 
magnetons fi* B — eh/2m* . 




FIG. 24: (a)-(e) contours of conditional densities \ip c (x,y)\ 2 
and phase of the conditional wavefunction 8 c (x,y) in gray- 
scale for L x = 5,L y = 10 four- minima QDM. See Fig. for 
details, (f) contours of total electron density of the three- 
vortex triplet state are plotted in the background and posi- 
tions of the vortices with the fixed electron in three different 
positions. 



anti-cross we see continuous increase of the magnetiza- 
tion before it starts to decrease again due to the contrac- 
tion of the electron density. It is interesting that after 
the bump the interacting magnetization has very similar 
dependence on the magnetic field as the non-interacting 
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FIG. 25: Density of the ground state at different magnetic 
field values for L x = 5, L y = 10 nm rectangular four-minima 
QDM. See Fig.[3for details. 



magnetization. At high enough magnetic field the elec- 
trons are localized into individual double dots and have 
single-particle properties and the spatial extents in the 
interacting and non-interacting systems are not very dif- 
ferent (see Eq. ©)• However, the electrons may move 
in a correlated way even if they are localized into one of 
the double dots. One should remember that there is also 
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FIG. 26: Most probable position in nm of singlet (S = 0) and 
triplet (S = 1) states for rectangular-symmetric four-minima 
(L x = 5, L y = 10 nm) QDM. Singlet ground state (magnetic 
field) region is marked with gray background color. 



the paramagnetic part in the magnetization, but this is 
constant for a given state, if the the angular momentum 
is a good quantum number, and does not depend on the 
magnetic field. Of course, in a non-circular symmetry 
the paramagnetic magnetization may not be constant as 
a function of magnetic field. 



D. Vortices of L x =5,L y = 10 nm QDM 

The phase information and the conditional densities of 
the rectangular four-minima QDM are shown in Fig. |2"H 
The most probable position lies now on the y axis. An- 
other possibility would be to have the most probable po- 
sition on a line connecting the two minima diagonally. 
However, as the other double dot is left with just one 
electron and the distance of L x = 5 nm is so small that 
the single-particle density is not localized to the minima 
of the double dot. So the most probable position is in 
the y axis. In rectangular symmetry correlations force 
the one electron to one of the double dots. With small 
L x s conditional density shows a peak at x = as in the 
non-interacting two-body density in Fig. |3| 

The conditional density becomes more localized as the 
magnetic field increases. The vortices appear sequen- 
tially in the QDM. The second and third singlet states 
(in Fig. |3|(c) and (e)) are not ground states as the sys- 
tem becomes spin-polarized after the first singlet-triplet 
transition. There is again a repulsion between the vor- 
tices. Interesting is to note that in (e) the two more 
distant vortices are positioned much closer to the fixed 
electron when compared to the double dot of Fig. 1121 (e). 
The white and dark regions near the borders of Fig. [31 
(e) show the shades of phase boundaries of more distant 
vortices (not visible in the figure). 

Vortex dynamics is studied in Fig. [21 (f ) for the three 
vortex triplet (at B = 6.5 T). The electron is fixed at 
three different positions and the total electron density of 
the same state is plotted in the background. One vortex, 
or Pauli vortex, is always on top of the fixed electron 



and the two additional vortices are symmetrically on the 
sides of the fixed electron. As the fixed electron is moved 
from the origin to the direction of the positive y axis, the 
vortices aside become closer to the fixed electron. 



E. Total electron density and the most probable 
positions of L x = 5, L y = 10 nm QDM 

Ground state electron densities in Fig. [31 show a lo- 
calization into two double dots as the magnetic field is 
increased. If the densities would be rotated by 90 de- 
grees they would resemble very much two-minima QDM 
(double dot) densities of Fig. [21 The smaller displace- 
ment (L x = 5 nm) in the four-minima QDM potential 
has a much smaller effect than the larger displacement 
(L y = 10 nm) because electrons localize into the double 
dots (with L x = 5 nm) separated from each other with 
a distance d — 2L y — 20 nm. Therefore electron density 
of rectangular four-minima QDM effectively resembles of 
that of a two- minima QDM (double dot). This is true 
for interacting two-electron system. 

The most probable positions of the lowest singlet and 
triplet states are shown in Fig. [31 Continuously chang- 
ing r* (i.e. no jumps) is due to anticrossing states where 
symmetry of a state (and also r*) changes continuously. 
Interesting is the strong suppression of the oscillations 
of r* at greater B. At large magnetic field the electrons 
localize into a distant double dots and interaction effects 
(like changing angular momentum states in parabolic 
QD) have a smaller impact on the properties of the two- 
electron system. The effect is quite different for parabolic 
QD and square-symmetric four-minima QDM where the 
localization of electron density is not so strong due to the 
nature of the confinement potential. 



VII. ANALYSIS OF THE RESULTS AND THEIR 
RELEVANCE TO EXPERIMENTS 

A. Role of symmetry in quantum dot confinement 

We start our analysis of the data presented above from 
the measurable quantities. One such observable is the 
total energy for ground and excited states as a function 
of the magnetic field, as well as the magnetization. To 
ease the comparison, these are collected in Fig. [31 The 
most striking feature is that the data of the the square- 
symmetric four-minima QDM is very similar to the one 
of the single parabolic dot. On the other hand, the data 
of the double dot resemble the one of a rectangular four- 
minima QDM. The reason behind the similarities of these 
pairs is the symmetry. The square-symmetric and circu- 
larly symmetric cases have higher symmetries than the 
rectangular ones. One can study this in detail by split- 
ting the total Hamiltonian of Eq. |QJ to two parts as 
H = Hq + Hi, where Ho is the Hamiltonian of the 
parabolic case (L x = L y = 0), and the impurity Hamilto- 
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FIG. 27: Three lowest singlet and triplet energy levels and 
magnetization for all studied quantum dot confinements. 



nian Hi contains the terms from finite L x and L y values, 
see Eq. J3J|. If we now have a high symmetry in the sys- 
tem, meaning L x — L yi one can see that the Hi does not 
couple the eigenstates of H that have a different sym- 
metry. On the other hand, in the case with L x ^ L y , Hi 
has a lower symmetry and more of the symmetries of Hq 
are broken. Due to this, states with different symmetry 
are coupled. This leads to anticrossings in the energies 
as seen in Fig. [23 where also the crossings of the high- 
symmetry cases are seen. One can estimate the strength 
of the symmetry-lowering part from the anti-crossing gap 
in energy, as in the point where the energies would cross, 
one has in the first approximation a Hamiltonian matrix: 



tt _ ( E E s \ 
\ E$ Eq J 



where Eq is the energy at middle of the gap, and 2E$ is 
the width of the gap. 

The magnetization curve for the low- and high- 
symmetry cases are also very different, see Fig. E7I 
A common feature in all these cases is the sharp in- 
crease in magnetization at the point where the total spin 
of the ground-state changes. On the other hand, the 
low-symmetry anticrossings of the energy result smooth 
changes in magnetization, whereas the high-symmetry 
data show sudden jumps. 

These findings indicate that it is rather difficult to ob- 
tain detailed information of the system based on the en- 
ergetics and the magnetization. The symmetry of the 
system can be extracted, but not much beyond that. For 
larger particle numbers, more and more of the correlation 
effects can be captured by the mean-field level. The effec- 
tive potential has, due to the Hartree potential, a higher 
symmetry than the mere external potential,- This results 
in less details to both energetics and the magnetization. 

A similar role of the symmetry can be seen on the non- 
measurable quantities, like the densities. On the other 
hand, the vortices are more delicate. This is because they 
depend linearly on the wave function, unlike densities and 
energies that are second order. 



B. Exchange of two spins 

The idea of double dot spin-swap operations in quan- 
tum computing™ lies in the coherent rotation of two ini- 
tially isolated spins. Starting with, say, spin-up elec- 
tron in the left and spin-down electron in the right dot: 
|T)i|l)fl an d rotating spins to opposite order: |J.)i|t)ii- 
These would be initial and final states of the system. Co- 
herent rotation between initial and final states requires 
entangled spin states like spin-singlet \S) = (\1)l\1)r — 
U>i|T)ji)/>/2and triplet |T > = (| T)i U) + U>i I T> i?)/ ^2 
whereas the other two triplet states |T+) = |T)l|T) R an d 
|T_) = \1)l\1)r are not conceivable as they have iden- 
tical spins. The states described above are only for the 
spin-part of the wave function but of course the spatial 
part, discussed extensively in this paper for singlet and 
triplet eigenstates, must be modified along with the ro- 
tation. 

In the simplified Heisenberg picture, rotation depends 
only on the singlet-triplet splitting energy, J = AE = 
E 1 * 1 * —E^, or exchange coupling of two spins. 14 Hubbard- 
type models can be used to study time evolution of a lit- 
tle bit more elaborate statesi 1 ^ 1 ^ To fully investigate the 
coherent rotation of two-electron system, it would be bet- 
ter to start with initially separated electrons (that can be 
constructed from many-body wave functions) and study 
the time evolution of the state in the exact many-body 
basis instead of using some simplified models. However, 
tuning J with dot-dot separation at small magnetic fields, 
the spin rotations can be quite safely modelled within the 
Heisenberg picture. At high magnetic fields, on the other 
hand, electrons in lateral double dots form finite quan- 
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turn Hall-like composite-particle states of electrons and 
flux quanta, as shown in this study and in previous stud- 
ies42il£ Therefore, the coherent spin rotations at high 
B may be quite different from zero B rotations, even if 
J could have exactly same value for high-_B and zevo-B 
states. These states are of course of great scientific in- 
terest as their own, but from the perspective of coherent 
two-spin rotations they may function quite differently as 
suggested by the Heisenberg or Hubbard models. Electric 
control of J with dot-dot separation may also be advan- 
tageous in other perspectives because magnetic fields are 
more difficult to apply locally. 

C. Comparisons to experiments 

In very recent experiments Petta et al. demonstrate 
a coherent rotation of two spins between singlet \S) and 
triplet \Tq) states in a lateral double dot device^ They 
start the operation from singlet state in single QD and 
then isolate two opposite spins in separated dots where 
the singlet-triplet splitting vanishes and no tunneling is 
allowed between the two dots. Coherent rotation is per- 
formed by bringing the two dots closer allowing small 
but finite energy splitting J between \S) and |Tq) states. 
Probability of finding singlet is measured as a function of 
gate operation time. Fig. (a) of this study shows cal- 
culated singlet-triplet splitting as a function of dot-dot 
separation and magnetic field. Following the zero mag- 
netic field line one and see how the singlet-triplet splitting 
decreases as a function of increasing dot-dot separation. 
Fixing B = 0.1 T, as in the experiment, for single dot 
(L = 0) we have J = 1.16 meV and for double dot at 
the greatest interdot distance studied (d = 2L = 40 nm) 
we have J = 0.16 meV. Petta et al~ find three times 
smaller value in single QD (J = 0.4 meV) as our calcula- 
tions. This is simply a consequence of different quantum 
dot confinement energy Hluq. 

Lee et al. studied experimentally singlet-triplet split- 
ting as a function of magnetic field in silicon two-electron 
double dot^ The measurements show very similar data 
as our results. The first singlet-triplet transition is re- 
solved clearly in the experiment in accordance to our 
calculations. Decreasing the coupling between the dots 
results in a small shift of J = to low fields as in our re- 
sults represented in Fig. |5J The high-field regime, where 
we would expect to find small J and even possibly a pos- 
itive J, which would correspond second singlet ground 



state, is not measured to very high field strengths. How- 
ever, to fully compare our results to the measurements on 
silicon double dots we should recalculate our data with 
silicon effective mass and dielectric constant. 

Brodsky et al. — were the first to measure ground state 
energy levels of a lateral two-electron double dot as a 
function of magnetic field. Even if they did not concen- 
trate on two-electron case particularly, the line for two 
electrons is clearly visible in their data showing also a 
kink indicating the singlet-triplet transition. 

Magnetization is very difficult to measure directly for 
just two electrons. However, indirect methods and direct 
methods with large arrays of individual few-electron dots 
may provide interesting experimental results^2£ Even if 
it is difficult to compare existing measurements of many 
electron dots to our two-electron system, the double dot 
measurements of Oosterkamp et al. ^ show similar type 
of anticrossings as our calculated magnetization curves. 



VIII. SUMMARY 

In summary, we have thoroughly studied different lat- 
eral two-electron quantum-dot molecules. All our exact 
diagonalization calculations were performed for closely 
coupled quantum dots. Many-body electron wave func- 
tions were allowed to extend over the whole system. We 
have analyzed how the physical properties change when 
a deviation or disorder is introduced in the confinement 
potential of the symmetric quantum dot. We have calcu- 
lated measurable quantities such as energy levels, singlet- 
triplet splitting, and magnetization as a function of mag- 
netic field strength. The measurable quantities were fur- 
ther analyzed by calculating non-measurable quantities 
such as phase vortices and conditional densities. We have 
also compared the properties of non-interacting electrons 
to interacting ones in quantum dot molecules to separate 
the effects of the non-circular confinement potential and 
interactions. 
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